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^ ■ Abstract 

We solve numerically the Boltzmann equation in the early universe in the pres- 
^ . ence of a constant electric field and find the electrical conductivity a in the range 

■ 1 MeV < T < Mw. The main contribution to a is shown to be due to leptonic 

interactions. For T < 100 MeV we find a ~ 0.76T while at T ~ Mw we obtain 
a ~ 6.7T 
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Recently, there has been much interest in primordial magnetic fields. Primordial 
fields might act as the seed field for the dynamo mechanism which may be responsible 
for the observed galactic magnetic fields |T]]. These have been measured both in the 
Milky Way and in other spiral galaxies, including their halos (which perhaps can 
be viewed as an indication of the primordial origin of the seed field). Typically the 
strength of the observed field is of the order of O{10^^) G. 

It has been suggested that primordial magnetic fields might be generated as a con- 
sequence of the early cosmic phase transitions, such as the electroweak or the QCD 
phase transitions 0. These would then be small scale random fields. They would be 
imprinted on the comoving plasma and dissipate very slowly [|| because the plasma 
of the early universe is known to have a very large electrical conductivity [Q. Indeed, 
on dimensional grounds alone one can argue that conductivity scales as a ~ T/a, 
where a is the coupling constant squared associated with the scattering processes of 
the charged particles in the plasma. As a consequence, the magnetic Reynolds num- 
ber of the universe, given by i? ~ vaH^^ where H = R/R ~ T^/Mpi is the Hubble 
parameter and v the bulk velocity, is huge. Despite its smallness, diffusion can play 
an important role in the subsequent evolution of the microscopic random magnetic 
field. The main issue here is just how the small scale field, the origin of which is mi- 
crophysics, gets amplified to a large scale seed field. In a recent paper Baym, Bodeker 
and McLerran |^ have suggested that in a first order electroweak phase transition, tur- 
bulence generated by the colliding shock fronts of the phase transition bubbles could 
give rise to large fields. More generally, a recent numerical simulation of the full 
magnetohydronamical equations in an expanding universe shows clearly a transfer of 
energy from small scales to large scales, albeit in 2 dimensions. Because of the very 
large Reynolds number, real numerical simulations are not presently realistic. How- 
ever, in 1+3 dimensions one may apply the so-called shell (or cascade) model, which 
has proved very succesful in studies of hydrodynamical turbulence, to the early uni- 
verse, with the result that the scale of the magnetic field fiuctuations indeed increases 
very rapidly @]. 

The previous estimates for the conductivity a have, however, been rough order 
of magnitude estimates only. The proper way to derive electrical conductivity of the 
early universe is to consider plasma in an electric field, and then use the Boltzmann 
equations to find the bulk motion, which serves to define conductivity through Ohm's 
law. To be able to do this, one first has to compute the collision integral, which is the 
weighted sum of the matrix elements of all the relevant processes leading to microscopic 
diffusion of the bulk motion. In the present paper we compute the collision integrals 
for leptons and hadrons (quarks) separately and solve the Boltzmann equations in 
the approximation where bulk velocity is much smaller than the thermal velocity of 
the plasma particles. This means that we treat the effect of the electric field as a 



perturbation in an otherwise homogeneous and isotropic background. 

The Boltzmann equation for the distribution function f{x{X),p{X)) of a charged 
particle in an electromagnetic field is found by requiring that along the world line, 
parametrized by A, the total change in / equals the collision integral. Using the 
equation of motion of a point particle with electric charge q, which in curved space 
reads 

% = -K,py' + qFj^p' , (1) 

one can cast the Boltzmann equation in the form 

f / + P • V/ - %,V^Rii - 2f A-f - E . p - ^ Cfe «)/ , (2, 

where C {p, t) is the collision integral. Here we have assumed a Robert son- Walker 
metric with a scale factor R for the background, and we have defined F° = —Ei and 
F/ = e^-'^-Bfc. We also prefer to use co- moving coordinates and define 

f{p,t) = J 5{po-{p'R' + m'y/')f{p,po,t)dpo . (3) 

Inserting this into Eq. (^) and integrating we find, making use of the local momentum 
defined as p = Rp, 

f-^f^-g^E^ = Cip,t). (4) 
dt R^ dp' ^dp' ^ ^ ^ 

Because of the assumed isotropy of the background, the terms p ■ V/ and E • pdf /dpo 
have been dropped. This is justified if one assumes, as we do, that the effect of the 
electromagnetic field on the distributions can be considered as a perturbation. From 
now on, we also drop the tildes for brevity. 

The collision integral encompasses all the scattering processes relevant for the dif- 
fusion of the charged test particle. For electromagnetism the leading processes are 
2 — > 2 reactions (we do not consider thermal effects such as the plasmon decay). Thus 

C(p, t) = - j dPbdP,dPdi2jc)^6\p + Pb-Pc- Pd)\M{Ab ^ cd)\^T, (5) 
Po -I 

where p is the four-momentum of the charged test particle A, and we have defined 
dPi = d^pi/ {{2TiY2E{pi)) . The factor JF contains the distribution functions of the 
initial and final states. For fermions one must include the Pauli blocking factors and 
for bosons the enhancement factors. Below are two examples of the jF-factor. The first 
is for the case when all the particles in the reaction are fermions, and in the second 
case particles h and c are bosons, such as is the case for Compton scattering: 

ni - f{p)][i - MPb)]fciPc)UPd) - [1 - /c(Pc)][i - UPd)]fip)MPb) , . 
1 [1 - f{p)][i + fMUMUPd) - [1 - + fc{Pc)]f{p)fM . ^ ' 



Above the time dependence is not shown exphcitly. 

Let us now assume that there exists a constant electric field in the early universe[] 
(constant in the sense that its coherence length is larger than the mean free path of 
charged particles) and treat its effect as a perturbation on the distribution of the test 
particle. We write f = fo + Sf, where /o is the equilibrium distribution and 6f the 
small perturbation. Inserting this into the collision integral Eq. (|^) results in 

C{p,t) = fdP,dPJP,{27Ty6\p + Pb-Pc-Pd) 

X \M\\[f,Jod[l- U + kb[l- kc][l- fm])5f ^C{p)6f{p,t). (7) 

Eq. (0) assumes that all the particles in the reaction are fermions; generalization 
to other cases is straightforward. 

We have assumed that all the particles b,c and d have an equilibrium distribution. 
Therefore, if also the test particle has an equilibrium distribution, the collision inte- 
gral vanishes. Thus only the perturbation term in Eq. (|^ survives. Note that the 
equilibrium distribution depends only on the momentum. 

Treating the electric field as a small perturbation, the Boltzmann equation can be 
linearized and reads 

with C{p) defined in Eq. (|^. Finally, assuming that ^ H\p\ = |p|, we may search 
for stationary perturbations, for which one easily obtains 

Sm = -J-^E^. (9) 

C{p) dp' 

The induced current density is given by 

j = 5fip)pd'p , (10) 

and conductivity a a, associated with a given particle species A, is defined by 

U = ^aE . (11) 



Thus the contribution of a single species A to conductivity is (Ta ~ 1/ \M{AX — > 
y)P, where the sum is over all the processes which scatter the test particle A. For the 
purpose of conductivity, we may view the mixture of different particle species, such as 

^Of course, we do not claim such a field actually exists but rather use it as a probe of the plasma 
properties. 



is found in the early universe, a multicomponent fluid. The flow of each component 
contributes to the total current and adds up to the total conductivity, which reads 

CTtot = '^o-A , (12) 
A 

where the sum is over all the relativistic charged species present in the thermal bath. 
Note that the total conductivity is dominated by the species that has the weakest 
interaction. This is because the weaker the interaction, the longer time the current 
flow is maintained. 

We consider the temperature interval 1 MeV < T < Mw, and make the simple, 
crude assumption that particles appear in the thermal bath only when temperature 
is greater than their mass. Thus below 100 MeV, for example, the only charged 
particles present are the electrons and positrons. When T > Tqcd, also the quarks 
should be counted in. Their main interactions are strong, so that their electromagnetic 
interactions may be neglected. The list of the relevant reactions for the leptons is 
presented in Table 1, and in Table 2 the same for the quarks, together with the matrix 
element squared (modulo factors of n and a). A technical point is that for the t- 
channel reactions there arises an infrared singularity, which is known to be regulated 
by thermal effects. These we approximate by giving fermions, photons and gluons 
a Debye mass in the t-channel and u-channel propagator; otherwise all fermions are 
assumed to be massless (consistency requires that the external particle masses are kept 
non-thermal). The thermal masses are given by 



2 ^ _j.2 _ 0.0115T2, 
' 8 

2 = - 0.251T^ 

6 

ml = 9!i3 + ^)T'^1.508{3 + ^)T\ (13) 

where Nf is the number of quark families present, the sums are over all particles with 
m < T, and we adopt the values = Airas — 1.508 and = inae — 0.0917, and 
Tqcd = 200 MeV. 



We have computed the perturbation Sf{p) numerically by evaluating the collision 
integral by a simple Monte Carlo simulation. The result agrees with the expected naive 
scaling law a ~ T. The form of Sf{p) is very similar for both leptons and quarks. As 
an example, in Fig. 1 we show dj /dp = qp^5f{p) for the electron at T < m^, where 



Table 1: Leptonic processes for e , and their matrix elements squared 
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Table 2: QCD-processes, and their matrix elements squared 



QQ 99 
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9s^ut 


(19 19 
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the peak at p ^ 3T is evident. In Fig. 1 we also compare Sf{p) with the equilibrium 
distribution fo{p) to demonstrate the region of validity of our calculation, where we 
have chosen E = lO^^T^ for definitess. For most part, 5/(p)//o is a constant, except 
when p — s> where Sf{p) vanishes. This behaviour is due to Compton scattering 
(and at higher temperatures also annihilation), for which the integrated matrix 
element squared gives rise to p-depencence slower than the one inherent in C{p). We 
should also point out that our result is sensitive to the the infrared cut-off provided 
by the thermal masses of the internal propagators, given by Eq. ([13D , and hence not 
strictly valid at very small p. This is because it is the t-channel reactions which 
give the dominant contribution to the collision integral Eq. (|^), and a straightforward 
numerical calculation gives an inversely linear dependence on the regulator in C{p)- 

With this caveat, the total leptonic conductivity, together with the total hadronic 
conductivity, is shown in Fig. 2 for the range 1 MeV < T < Mw, above which one 
would have to account for the charged W's. However, purely electromagnetic pro- 
cesses continue to dominate conductivity also well beyond the Z-pole. The reader may 
compare the result presented in Fig. 2 with the text-book estimate Q for relativistic 
electron gas scattering off heavy ions, which yields a ~ 14T. The steps in the leptonic 
a reflect both the appearance of new leptons in the thermal bath, as well as change in 
the regulator (thermal photon mass) due to the appearance of new quarks and leptons. 

Introducing a thermal mass into the propagator is only an approximation of the 
full thermal screening, which includes also initial state and vertex corrections. Such 
corrections should remove the spurious infrared singularities in C(j)). Our approxima- 




Figure 1: dj/dp and 5///o for the electron at T < m^. Here E = 10 ^T^. 
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Figure 2: cr/T as a function of temperature. 



tion should, however, be sufficient to produce a reliable order-of-magnitude estimate 
of the electrical conductivity in the very early universe. 
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